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Abstract 


We have searched for the rare decay of the eta meson r] e~^e~ using the 
CLEO II detector. The r/’s were produced in e'’"e“ collisions with 10 GeV 
center-of-mass energy at the Cornell Electron Storage Ring (CESR). We find 
with 90% confidence the upper limit on the branching fraction R (?7 —> e'’'e“) < 
7.7 X 10“^. The application of conventional elementary particle theory to this 
decay predicts a branching fraction of about 10“®. 
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I. INTRODUCTION 


We have used the CLEO II detector at the Cornell Electron Storage Ring (CESR) to 
study about 2 x 10 ^ events of the form e’''e“ —>■ hadrons to search for the rare decay mode 
rj ^ e~^ + e~. We have not found this decay, only a 90% conhdence upper limit of B{ri ^ 
e~^e~) < 7.7 x 10“^. As discussed in Sec. II on the conventional theory of this decay, the 
predicted branching fraction is about 10“®. An observation of a signal above this level could 
be evidence for an unconventional process which enhances the rj —>■ e’''e“ decay rate. 

The plan of this paper is as follows. Section II gives the conventional theory for the decay 
and the predicted relationship between i?(p —»• and the measured branching fraction 

R(p —Then Sec. Ill describes the data and the calculation of the number of r^’s 
produced in the 2 x 10^ hadron events. Section IV describes the method used to search for 
r] —> e+e“ decays. Finally Sec. V contains the calculation of the upper limit on i?(p —> e’''e“), 
a discussion of the errors, and some general remarks on the search. 


II. CONVENTIONAL THEORY FOR r] e+e" AND t] p+p” 

The rj is massive enough (547 MeV/c^) to decay via rj p+p“ as well as via rj —>■ 

Figure shows the decay mechanism according to conventional theory Si- The decay 
matrix element for rj into two virtual photons, represented by the cross hatched circle, is 
difficult to calculate from Erst principles. Indeed the same problem occurs in the study of 
the decay 7 r° —>• e+e“ l]!],^. 
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FIG. 1. Feynman diagram for the conventional theory of the decays rj —> e^e and rj 
I represents an e or a p. 




Landsberg p| has reviewed the conventional theory for the decay P Pi where P is 
a pseudoscalar meson, and he gives the formula for the branching fraction 


B{P l+r) = B{P 77)2aVs[|X|2 + |y|2]. 


( 1 ) 


Here B{P 77 ) is the branching fraction for P —>■ 77 , a is the hne structure constant, 
r = rrii/mp where mi and mp are the I and P masses, and s = (1 — 4 r^) 2 . Y is proportional 
to the imaginary part of the decay amplitude and has the explicit form 


|V| = lln 

s \ 2r 


( 2 ) 


4 







X is proportional to the real part of the decay amplitude, and is difficult to calculate with 
certainty. Of course even if X is calculated precisely, Eq. o is still not a basic formula 
because B{P ^ 77 ) must be taken from experiment. Nevertheless Eq. (|^ with X = 0 gives 
a lower limit on B{P ^ namely 


B{P /+/-) > B{P 77 ) 


2a^r 


2^2 


In 


1 + s 
2 r 


(3) 


Using the rj —>■ branching fraction from the Particle Data Group |p , 


B{v 77 ) = 0.393 ± 0.003 


(4) 


and Eq. (j^) we can calculate 

B{t] ll^ll~)min = 4.4 X 10“®. 


(5) 


The measured branching fraction |^|-|^, 

£?('/; ^ 1-^'*'li~)measur^d = (5.7 ± 0.8) X lO"*" 


( 6 ) 


is consistent with this limit. If we assume that the ratio |Xp/|Up is the same for rj —>■ e'''e 
and 7] —>• then 


B{ri —>• e’''e ) 

B{v /iV“) 



(7) 


where the e and p subscripts refer to me and in r = mi/mp. The dominant term, 
{j'el'fiiY = {me/mfj,Y, is due to helicity suppression. From Eq. (|^), 


B{r] 


e^e 


= 4.05 X 10 


-4 


B{v ) 

Eq. (|) also holds for Bij} ^ e^e~)min/B{rj Using this and Eq. 


( 8 ) 


Biv 


e e 


= 1.8 X 10 


-9 


(9) 


Finally, combining Eqs. (H) and (§), 

B{7] e+e“) 2.3 x 10"®. (10) 

This estimated branching fraction is based on the assumption that |Xp/|yp is the same 
for the e+e“ and p+zi” decays. An unknown process present in the decay 7 —> e+e“, but 
not in the decay 7 —>• could result in a value of |Xp/|yp much larger in the e+e“ 

mode, and thus a signal larger than the above limit. 
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III. DATA AND NUMBER OF t^’S PRODUCED 


A. Detector and Data 

We used data collected by the CLEO II detector p at CESR. The components of the 
detector which are most critical to this study are the three concentric cylindrical drift cham¬ 
bers occupying the space 4 cm to 95 cm radially from the beam axis, comprising a 67- 
layer charged-particle tracking system which is immersed in a 1.5 Tesla solenoidal magnetic 
field. The momentum p, in GeV, of charged particles is measured with a resolution of 
<Jp/p(%) ~ [(0.15p)^ -|- (0.5)^] 2 0.5% for the electrons we see in rj decay. In addition, 

ionization loss {dE/dx) is measured in the 51-layer main drift chamber with a resolution of 
6-7%. 

Also important for the measurements reported here is an electromagnetic calorimeter 
consisting of 7800 thallium-doped Csl crystals. These crystals, each of dimension ~ 5 cm x 
5 cm X 30 cm, surround the tracking volume, convering 98% of the full solid angle. Forming 
the barrel region of the calorimeter, 6144 tapered crystals are arrayed just inside the magnet 
coil at a radius of ~ 1 m in a projective cylindrical geometry, covering 82% of the solid angle. 
The remaining crystals are rectangular, and are oriented axially in two end caps, overlapping 
in solid angle with the ends of the barrel. The barrel region of the calorimeter achieves 
energy and angular resolutions for electromagnetically showering particles of aE/E{%) = 
0.35/i7°''^^-l-1.9 —O.IE and cT</,(mrad) = 2.8/\/E + 2.5 {E in GeV), respectively. The resulting 
photon energy and direction information provided by this system is used to reconstruct 
p —77 decays to determine our rj sample size as described below. The crystals are also 
used in distinguishing electrons from pions. When combined with the tracking and dE/dx 
information from the drift chambers, the misidentihcation of pions is limited to less than 
one percent. 

We used 3.11 fb“^ of data at the T(4S) resonance, 10.57 GeV, and 1.69 fb“^ of data 
below the T(4S) resonance at 10.53 GeV. The samples of on resonance and off resonance 
hadronic events were: 


Nhadion) = 1.78 X 10^ (11) 

7V;,,rf(off) = 0.797 X 10". (12) 

We note that 75% of the on resonance events are of the same type as the off resonance 
(also called continuum) events. The remaining 25% of the on resonance events are from 
e’''e“ —BB production. We found that a major source of contamination in the search for 
rj —> e’''e“ decay is one real electron from the semileptonic decay of the B or B plus one 
false electron from a pion from the B or B decay chain. We substantially reduce the number 
of BB events by applying the following selection criterion to the on resonance events. We 
require R 2 > 0.3, where R 2 = H 2 /HQ and Ho,H 2 are the zeroth and second Fox-Wolfram 
moments, respectively . This selectively removes almost all of the more spherically-shaped 
events (see Appendix A for clarihcation). Thus when discussing our total sample size, we 
emphasize the number of continuum events in our data sample: 

Nhadion) = 1.34 X 10" (13) 
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NhadioS) = 0.797 X 10 ". 


( 14 ) 


Note that this requirement also removes about hfty percent of the more collimated continuum 
events; this effect is included in the efficiencies we determine below. 

We determined the number of tj mesons in the events in Eqs. (|^) and ( p!4[) using the 
following procedure. We selected pairs of photons from the decay 


r] 77 . (15) 

We then used continuum e’''e“ —> hadron events generated from Monte Carlo and our 
simulation of the properties of the CLEO II detector to determine the efficiency for detecting 
V 77- We also calculated the average number of r/’s produced per event. We now give the 
details of the procedure. 


B. Observed 7 — > 77 Decays 

In developing the criteria for selecting photon pairs from 7 —>■ 77 we kept in mind that we 
would be looking for 7 —> e+e“. As much as possible we chose the same selection criteria for 
7 pairs as we would use for the e+e“ pairs. In this way any uncertainties in the Monte Carlo 
modeling of the events will be applied to both decays and thus cancel out of the analysis. 
For example, since a crucial identihcation signal for both 7 ’s and e’s is an electromagnetic 
shower in the calorimeter, and since a minimum shower energy of 0.4 GeV is required for 
precise e — vr separation, we set the minimum shower energy at 0.4 GeV for 7 ’s as well. 

Only e’''e“ —>■ hadron events were used by requiring the following selection criteria. First, 
the events must have at least hve charged tracks each of momentum greater than 225 MeV/c. 
Second, non-annihilation events such as those from beam-gas or beam-wall interactions were 
rejected. Third, for events taken at beam energies corresponding to the T(4S') resonance, 
BB events were removed by the R 2 cut described previously. Further selection criteria for 
identifying 7 ’s from 7 —>■ 77 were: 

a. The electromagnetic shower must occur in the main (barrel) portion of the electromag¬ 
netic calorimeter. Specihcally, jcos^l < 0.71 where 6 is the angle between the 7 direction 
and the beam axis. 

b. The angle between the shower and the nearest charged particle track must be larger 
than 20°. 

c. The shower energy must be larger than 0.4 GeV. 

d. The pattern of energy deposition of the shower in the crystals must be characteristic 
of a single photon (two or more clusters of deposited energy indicate randomly overlapping 
photons or a high-momentum 7r°). 

e. The shower must not appear to be a fragment from another shower or from a charged 
pion interaction. 

f. If two 7 ’s have an invariant mass within 2.5 ct (12.5 MeV/c^) of the 7 r° mass, both 7 ’s 
were discarded. 

Using 7 ’s selected with these six criteria we then considered all combinations of pairs 
subject to the condition that 


P 1 +P 2 I > 0.8 GeV/c 


( 16 ) 
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where pi and p 2 are the vector momenta of the two 7 ’s. We set this condition to rednce the 
nnmber of random pairings of 7 ’s. 



FIG. 2. Measured 7 pair invariant mass distribution for (a) on resonance events and (b) off 
resonance events. 

Figure shows the invariant mass spectra for 7 pairs for the on and off resonance data. 
The 7 mass peaks at 546 MeV/c^ have the widths expected from the properties of the CLEO 
II detector, namely a = 12.8 MeV/c^. Fits to these spectra gave the following number of 
observed 7 —>• 77 decays: 


iV^^.^.^(on) = (7.59 ± 0.12) x 10^ (17) 

= (6.57 ± 0.06) x lOh (18) 

Hence we observed about 1.42 x 10® 7 —> 77 decays. 


C. Efficiency for Detecting 7 —> 77 

In determining the efficiency for detecting 7 —77 events we restricted ourselved to 
simulated e’''e“ —> hadron events in the continuum. Using the selection criteria described in 
Secs. IIIA and IIIB and the known number of 7 ’s in our Monte Carlo sample, we calculated 










the efficiency for on resonance (with the R 2 cut) and off resonance (without the R 2 cut). 
Since from Eqs. (]^) and (M) we see that 37.5% of our events are off resonance and 62.5% 
of our events are on resonance, we use these fractions to hnd a weighted mean efficiency for 
V 77, 


= (4.80 ± 0.05)% (19) 

where the error is statistical. The systematic error, which is substantially larger, will be 
discussed in Sec. V. We also found good agreement between the observed 77 —>• 77 events 
and the simulated 7 ^ 77 events with respect to the 7 momentum spectrum, the rj an¬ 
gular distribution, and the 7 pair invariant mass distribution of both the rj peak and the 
background. 

One major reason for the small efficiency in Eq. ([I^) is that the 7 momentum, is 
required to be larger than 0.8 GeV/c, but most p’s are produced at smaller momenta. 
Also responsible for the small efficiency is the requirement that the 7 shower energy in the 
calorimeter, E^, be greater than 0.4 GeV. Smaller lower limits on Pn and would increase 
substantially, but would result in large increases in background for rj —>• 


IV. SEARCH FOR p ^ e+e" DECAYS 


As already noted we limited our systematic uncertainties, particularly our dependence on 
simulated event sets, by using as much as possible the same selection criteria for p ^ e~^e~ 
events as we used for the observed p —>• 77 events. Secs. IIIB and IIIG. 

Beginning with the same event sample in Eqs. (|T^ and (|T^ we again require hve charged 
tracks, classihcation as an annihilation event, and R 2 > 0.3 for the on resonance events. We 
then looked for showers associated with a charged track which met the following criteria: 

a. The shower angle 6 must satisfy | cos^l < 0.71. 

b. The shower energy must be larger than 0.4 GeV. 

c. The energy deposition of the shower in the crystals must be characteristic of a single 
photon. 

d. The shower must not appear to be a fragment from another shower or from a charged 
pion interaction. 

e. We require IP 1 +P 2 I > 0.8 GeV/c. 

The charged particle track had to satisfy our standard criteria for a track from the 
primary interaction vertex, namely: 

f. The track had to be of good quality, as identihed by the GLEO software tracking 
algorithms. 

g. The distance of closest approach of the track to the beam line had to be less than 5 
mm. 


h. The distance of closest approach of the track to the event vertex measured parallel to 
the beam line had to be less than 50 mm. 

Next the track had to be identihed as an electron using a standard GLEO algorithm 
combining E/p, shower shape, and several other parameters. The algorithm has an efficiency 
of greater than 90%, with the exact efficiency depending on Pe, and a fake rate from charged 
pions of about 0.5%. 
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FIG. 3. Invariant mass spectrum of e"*'e“ pairs found in the search for r] e"*“e“ decays. There 
is no peak at the rj mass of 0.547 GeV/c^. The dashed curve shows the 90% upper limit for the 
r] —> e^e“ signal plus background. 

We then calculated the invariant mass of every e’''e“ combination in the events of Eqs. 
0) and (H). Figure ^ shows the spectrum in the mass range of 0.5 to 0.6 GeV/c^. There 
is no peak at the rj mass of 547 MeV/c^. A study of simulated rj —> decays showed 

that a peak would have a a of about 5 MeV/c^. 

V. CALCULATION OF UPPER LIMIT, ERRORS, AND FINAL REMARKS 
A. Calculation of Upper Limit on B{r] ^ e+e“) and Errors 

To determine the upper limit on the branching fraction B{'q ^ we have to know 

the efficiency for the detection of ?7 —*• e+e“ decays using the criteria in Sec. IV. We generated 
simulated 7] —*• e+e“ decays, applied these criteria and found the total efficiency by taking a 
weighted mean of on resonance and off resonance data as in Sec. IIIC. The mean efficiency 
was found to be 


= (5.22 ± 0.31)% 


( 20 ) 


where the error is statistical. 
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We used the same simulated events to parameterize the shape of a hypothetical rj —>• e+e“ 
peak, resulting in a mass of (545.9 ± 0.1) MeV/c^ and a width (a) of (5.4 ± 0.1) MeV/c^, 
where the errors are statistical. This shape was then used to £t the data of Fig. to hnd a 
90% confidence upper limit on the number of p —>■ e+e“ decays, iV^^e+e-- Varying the mean 
and a of the £t by one standard deviation, and alternately applying linear and quadratic 
background functions, gave a range of values for V^^e+e- from 18.4 to 27.1 events. We used 
the most conservative of these hts and concluded 

< 27.1. (21) 

It is at this point that we must consider the major sources of systematic error in our 
analysis. As we have already minimized the uncertainty from the electromagnetic calorimeter 
selection criteria, our remaining sources of systematic error stem from our tracking, particle 
idenhcation, and photon detection efficiencies, as well as the uncertainties in 
^'rj^e+e-^ Bij} ^ 77 ). The total systematic uncertainty of 10.7% is calculated in Table 
|. In our hnal calculation we reduce our mean efficiency by this amount, yielding: 

£^^e+e- = 4.66%. (22) 

Finally, we normalize to the branching fraction of p — 77 , 

B{7] —>• e’''e“) = X X B{ri 77 ) (23) 

I^rf —^77 e~ 

and arrive at our 90% confidence upper limit, 

B{r] e+e") < 7.7 x 10“^ (24) 

This upper limit is indicated by the dashed curve in Fig. ^ 


TABLE I. Summary of Systematic Uncertainty 


Source 

Uncertainty 

Tracking efficiency 

1 % per e candidate 

Electron ID efficiency 

3% per electron 

Photon detection efficiency 

3% per photon 

(stat.) 

1.4% 

—^77 (stat.) 

1 .0% 

(stat.) 

6 .0% 

B{v 77) 

0.7% 

Total 

10.7% 
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B. Final Remarks 


Our limit of B{ri ^ e~^e~) < 7.7 x 10“^ agrees with and improves upon the upper limit of 
2 X 10“^ found by White et al. 0; both of the conhdence levels are 90%. White et al. used 
the reaction p + d —^^He +77 to produce r/’s combined with a two-arm counter telescope to 
search for the rj —>■ e+e“ decay. In the course of concluding our analysis we have considered 
if improvements could be made in our method. 

We note from Eqs. ( 0 ) and ( p! 8 D that about 1.4 x 10^ r; —>• 77 events were observed. 
Since £,(^77 and s^^e+e- are about the same, we should have been able to investigate a 
B{t] ^ e'^e~) of the order 

B{r] e+e") ~ ^ 77 ) ~ 6 x 10“® (25) 

if there were no background events. In the future when the number of detected e’''e“ —^ 
hadron events increase twofold or more at CESR, and at the B-factories now under con¬ 
struction, one might hope to achieve a sensitivity of 10 “®. 

However there is a background primarily from pairs containing one true electron and 
one pion misidentihed as an electron. Removing BB decays aided us somewhat, but unless 
this background is further reduced, sensitivities of 10“® to 10“® for B{r] ^ e~^e~) cannot be 
achieved by our method. The RICH detector that will be installed for CLEO HI should help 
improve the pion-electron separation. In addition, substantial improvement in sensitivity will 
probably be achieved using hxed target 7 production via hadronic collisions and specially 
designed electronic detectors. 
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APPENDIX A: VALIDITY OF SUPRESSING BB 

Our efficiencies for the on resonance and off resonance searches for r] 77 , i.e. with and 
without the R 2 cut, respectively, are 

e^^^^ion) = (3.75 ± 0.05)% (Al) 


£^^.^.^(off) = (6.57 ±0.07)%. 


(A2) 


Combining Eq. (^) with the number of events in our off resonance sample [Eq. ([l^) ], the 
number of 7 ’s observed [Eq. (P!B|)], and B{t] ^ ll)i we hnd the number of r/’s produced per 
continuum event to be: 


u^(cont.) = 0.315 ± 0.006. (A3) 

Multiplying this result by the number of on resonance continuum events [Eq. (|^], the on 


resonance efficiency [Eq. (^)], and B{rj ^ 77 ) gives the number of rj 
should expect to see from on resonance continuum events. 


77 decays we 


A'^^.^.^(on) = (6.50 ± 0.25) x lOh (A4) 

Subtracting this from the number of decays we do see, Eq. [(0)], this leaves about 10900 ± 
2800 7 —77 decays that must come from BB events, or about (8±2)% of the total number 
of on and off resonance decays we observe. Thus it seems that our assumption that all BB 
events are supressed leads to an small overestimate of 

However, one must remember that and Aj^^e+e- will have almost equal propor¬ 

tional contributions from BB events. In fact the only discrepancy between the two con¬ 
tributions will be due to slightly different acceptances, caused by the differing angular and 
momentum distributions of the two channels. Since in calculating our hnal limit we take 
the ratio the BB contributions will almost entirely cancel. The remaining 

effect will be much smaller than 8 %, and thus negligible for an upper limit. 
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